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The present paper deals with the equihbrium properties of spin- 2 non- 
uniform anisotropic XY chain in transverse field that is defined by the 
Hamiltonian 

N N-l 

^js^j + E (jr^js^i + (1) 
j=i j=i 

The quantities of main interest are the time-dependent correlation func- 
tions < s'^{t)s'^^^ >, < (...) Sp[e-^-f^(...)]/Spe-^-f^, sj{t) = e'^'s'^e-'^^ 
and dynamical susceptibilities Xap{^,uj) = Ene^^'^/o^ c^te^('^+")*i < [s'^{t), 

The system considered in partial cases has well-known ferroelectric inter- 
pretation in connection with some quasi-one-dimensional hydrogen-bonded 
ferroelectric crystals like Cs{Hi^xDx)2P0i, PbHi_xDxP04 etc. Indeed, 
in such compounds the hydrogen bonds, that play an important role for 
the ferroelectricity, run along chains, the intrachain interactions are much 
stronger than the interchain ones, and a suitable model for describing pro- 
tons behaviour in a chain is spin-i Ising model in transverse field, which 
is contained in (1) (Jj^ = 0, Jj^ is interaction between neighbouring hy- 
drogen bonds along chain and flj describes proton tunneling motion). The 
non-uniformity may be caused physically by deuterization of some hydro- 
gen bonds. Under simplest assumption one believes that the substitution 
H ^ D changes significantly the tunneling probability and practically does 
not affect the intcrsitc interaction parameters. For this case one should put 
jxx _ J g^j^j suppose that flj varies randomly from site to site taking value 
ft with probability 1 — x and zero value with probability x, that is, one 
should consider random chain (1) with the probability distribution density 

N 

n„ jf , jf , ...) = nM(^.) + (1 - - ^M^r - mJT)- (2) 



Statistical mechanics calculations for this model bear strong relevance for 
paraelectric phase properties of mentioned crystals and have been performed 
mainly in the limiting cases x — 1 and x — (see e.g.^~^ and references 
therein). Such studies were encouraged by the fact that the model (1) with 
the help of Jordan- Wigner transformation can be rewritten as a system of 
non- interacting fermions and thus a lot of result can be obtained exactly^. 
Nevertheless, there are serious difficulties in calculation of longitudinal time- 
dependent correlation function < s^{t)sj_^^ > and therefore the dynamical 
susceptibility Xxxi^,'-^), and the complete solution even for perfect (non- 
random) case has never been found. The results obtained within different 
approximate approaches have serious shortcomings. 

In our previous papers^"^ we suggested a new numerical approach for 
evaluation of thermodynamical properties and both equal- time and different- 
time spin correlation functions for the model defined by (1). In the present 
paper the proposed approach is described in more detail and is applied for 
examining the influence of disorder on the average transverse dynamical 
susceptibility of spin-| Ising chain in random transverse field (1), (2). Such 
results, up to the best authors' knowledge, are obtained for the first time. 

In order to obtain the equilibrium statistical properties of the system 
in question (1) one should first perform Jordan- Wigner transformation from 
operators sf = sj±is| to Fermi operators Ci = , Cj = nn=i(~2s^)sj, j — 
2, ...,N and then to diagonalize the obtained quadratic form by a canonical 
transformation 77^ = J2jLi{gkjCj+hkjCj') with the outcome H — J2k=i ^kiVk^k — 
^). The elementary excitation spectrum and the coefficients gkn = 



($fen + ^fcn)/2, hkn = {^kn - '^kn)/'^ are determined from 

N 

J=l 

TV 

Afe^'fen = E (A - B) jn, (3) 

jxx_^jyy jxx j^jyy jxx_jyy 

where Aij = fliSij + ' 4 ' (5^,^+1 + '"'4 -B^j = ' 4 ' - 

jpa; —jyy 

'"^4 Eqs. (3) reduce to standard problems 

TV 

A^,$,„ = $fe,-[(A - B)(A + B)],„, 

AT 

A^^fen = E + B)(A - B)],„ (4) 

for N X N five diagonal banded matrices (A ^ B)(A ± B). 

Next it is necessary to express the quantities of interest through A^, 
^kj, '^kj- Since the partition function Z{P,N) = IlfcLiScosh^ involves 
all Af,, the knowledge of eigenvalues of matrices (A =]= B)(A ± B) permits 
one to obtain the thermodynamical properties of the model in question 
(1). Reformulating spin operators in terms of operators = + Cj — 

EZ=l^mj(Vm + Vm) and ipj = C+ - Cj = Em=l*mj(^m " Vm) , that iS, 



5j = 4 = UlJi{v>;M/(^')' = -^/^7/2, noting 



that the calculation of spin correlation functions reduces to exploiting Wick- 
Bloch-de Dominicis theorem, and calculating the elementary contractions 

< ^tit)^m E*p^-*p- — Aa, ^ 

p=i cosh ^ 



sinh(-iApt + ^) 



^ cosh 



^ coshfiA t — ^) 

< <pj{t)v;n >= - E — ' ^ (5) 

p=i cosh ^ 

one finds that ^kp, "^kg for all are involved in final expressions for equal- 
time and different-time spin correlation functions. Further numerical com- 
putation may be performed in a following way: first, to solve the eigenvalues 
and eigenvectors problem for matrix (A -|- B)(A — B) (4) obtaining in result 
A^ and ^^j, and second, having "^kj and A^ = ^/^, to find ^kj from Eq.(3); 
thermodynamics and spin correlation functions are expressed via the sought 
quantities. It should be stressed that here appears only the eigenvalues and 
eigenvectors problem for N x N five diagonal banded matrix that is the 
remarkable peculiarity of the system under consideration (1). This is a key 
difference of our approach in comparison with other finite-chain calculations 
(reported e.g. in^) where the numerical diagonalization of the Hamiltonian 
(1) in the Hilbert space of dimension 2^ is performed. 

The remainder of the paper deals with numerical calculation of trans- 
verse dynamical susceptibility Xzz{^,^) for spin-i Ising chain in random 
transverse field (1), (2). Acting as described above one finds that 4 < 

Sji^)Sj+n > = < Vjit)vj{t)ipf+^ip-+^ > = < ip+it)ip-{t) >< vj+nVj+n > 
- < "fti^W^+n >< Vjit)Vj+n > + < ^jitWj+n >< ^jit^+n >, < 

Sj_^_nSj(t) >=< Sj^^{—t)sj >. Computing then numerically the elemen- 
tary contractions (5) involved in the correlation functions < •Sioo(^)'5ioo+n 
< •5ioo+n( "0*100 > for the random chain of 200 spins with J — —2, Q — 1 
at (3 — 20, performing numerically the integration over t and summation 
over n, and averaging over the realization (typically few hundred) we have 
obtained the results some of which are presented in Fig.l. As the concen- 
tration of the sites with transverse field fl increases Xzzi^,^) rebuilds from 



the generally known Ising-like type behaviour to the behaviour inherent to 
Ising model in transverse field derived for the first time with periodic bound- 
ary conditions imposed in^. Note, that the depicted frequency shapes lose 
their symmetry and develop into smooth curves ai x — 0. The obtained 
susceptibilities exhibit a lot of structure. The detailed structure is induced 
by the disordered arrangement of two values of transverse field and Q = 1. 
Considering at first the case of small concentrations of one can find that 
each well-defined peak is associated with the susceptibility of certain chain 
that is determined by values of transverse field only in the local environment 
of spin at site j = 100 (see Figs. 2-4). While x decreases the number of possi- 
ble local structures (and thus the number of peaks) increases and the peaks 
appear almost at all frequencies. Nevertheless, since the difference in their 
magnitudes conditioned by the probability of their appearance is large, even 
at rather small x it is still possible to recognize the peaks corresponding to 
some simple local structures. Completely smooth curves appear only in the 
limiting case x = 0. It should be added that the similar calculations for 
some other types of disorder yield smooth frequency shapes. Finally, one 
should remind that the obtained results are exact only for finite chains and 
it is necessary to be careful while making extrapolations to A?" = oo, but, 
apparently, the picture qualitatively will remain the same. 

In conclusion it should be underlined, that the suggested method appears 
to be of great use in study of dynamical properties and random versions of 
spin- 1 XY chains and the derived results may be useful for interpretation 
of observable data in dynamical experiments on quasi-one-dimensional com- 
pounds that can be described in frames of model (1). 

The present paper was presented at Plenary Session at Ukrainian-Polish 



& East-European Workshop on Ferroclcctricity and Phase Transitions (Septem- 
ber 18-24, 1994, Uzhgorod - V.Remety, Ukraine)^0. The authors would hke 
to express their gratitude to participants of the Workshop for interest in this 
paper and discussions. 
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Figure captions 

Fig.l. Frequency-dependent real and imaginary parts of transverse suscep- 
tibility Xzzlas, f^) for spin-i Ising chain in random transverse field (1), (2) at 
ae = for different values of concentration x; e — 0.01. 

Fig. 2. Xzzi^i^) for certain spin-^ Ising chains. It is determined by values of 
transverse field in local environment of spin at site j — 100 (no transverse 
fields) and not denoted Q^s do not influence it; e = 0.01. 

Fig. 3. The same as in Fig. 2 for one transverse field in local environment of 
spin at site j — 100; e = 0.01. 

Fig.4. The same as in Fig. 2 for two transverse fields in local environment of 
spin at site j — 100; e = 0.01. 
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